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In  this  lecture  I  am  going  to  outline  some  of  the  early  ideas  about 
transonic  flow  with  which  Hans  U.  Liepmann  was  associated  and  show  how  they 
look  in  light  of  later  developments. 

The  earliest  ideas  on  transonic  flow  can  be  traced  back  to  studies  on 
gasdynamics.  Chaplygin's  paper*  in  1902  treated  planar  gas  jets  by  the 
hodograph  method.  Modern  theoretical  work  connected  to  aeronautics  dates  from 
the  papers  of  von  Karman^,  Guderley^,  and  Frankl  who  all  derived  the  approx¬ 
imate  equation  of  transonic  flow  around  1946.  Early  experimental  work  was 


carried  out  by  Stack  and  Dryden  at  NASA  Langley  in  the  early  '40’s.  Karnan's 

paper  represented  the  velocity  potential  ♦  for  flow  past  an  airfoil  (as  in 

Fig.  1)  as  a  uniform  flow  at  the  critical  speed  a*  plus  a  small  disturbance 
4»(x,y) 

«  =  a*  x  +  $(x,y)  . 

Assuming  that,  because  of  the  transonic  nature  of  the  flow,  JL  »  JL 
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Karman  derived  the  equation 
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This  basic  nonlinear  equation  of  changing  type  is  at  the  heart  of  all 
transonic  theory.  Karman  also  noted  the  similarity  parameter 
(effectively) 

l-M* 

K  * _  5  *  t/c  ,  M  «  iL  *  Mach  number  at  infinity  *  flight  speed  (2) 

6 2/ 3  •a.  sound  speed 

and  gave  scaling  rules  for  airfoil  flows.  Liepmann,  Ashkenas,  and  Cole*  gave 
a  more  detailed  derivation  as  well  showing  how  including  compression  viscosity 
can  yield  a  viscous  transonic  equation  which  can  have  smooth  solution  even  when 
shocks  are  present.  The  shocks  were  shown  to  be  thin,  Res  *  Reynolds'  number  of 
the  shock  ■  1. 

Topics  of  considerable  interest  were  the  possibility  of  obtaining  shock 
free  mixed  subsonic-supersonic  flows  by  the  use  of  the  hodograph  method  and  the 
physical  significance  of  the  limiting  line;  also  the  effects  of  the  viscous 
boundary  layer  on  the  inviscid  flow,  especially  in  shock  wave-boundary  layer 
interactions.  These  are  discussed  below. 

Another  significant  achievement  of  von  Karman's  was  to  bring  Hans  W. 

Liepmann  to  Caltech  where  he  stimulated  and  influenced  a  whole  generation  of 
students  and  created  an  early  interest  in  transonic  flow. 

1 1 .  Transonic  Small -Disturbance  Theory 

The  equations  of  transonic  small -disturbance  (TSD)  theory  can  be  regarded 
as  part  of  systematic  limit  process  expansion.  The  starting  point  is  the  Euler 
equation  of  inviscid  compressible  flow  and  the  Rankine-Hugoniot  jump  conditions 
for  shock  waves,  including  the  condition  that  the  entropy  increase.  A  typical 
geometry  is  shown  (in  transonic  coordinates)  in  Fig.  2  where  a  vortex  sheet 
trails  downstream  of  the  lifting  wing.  The  asymptotic  expansion  has  the  formS 
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^  »  (1  +  6 2/ 3  u(x,y,z;K)  +  5^/3  U£  ♦  ***)  +  S  v  +  6^/3  V£ 

♦  V  =  (v,w)  (3) 

p  ^  ^  ^  ^ 2/ ^  p  ^  iii  0  ^  ^  ^  ^2/3  ^  iii 

The  limit  process  has  4  ♦  0,  M?  »  1  -  tCS^/ 3  ♦  \  with  (x,y  *  ftl/3y,  z  *  fil/^z, 

« 

1-M2 

K  *  _ -  fixed).  Lengths  are  measured  in  terms  of  the  typical  wing  chord  c. 

6  2/  3 

The  span  b  should  grow  as  6  ■*■  0  such  that  B  *  b6*/3  js  fixed.  Timman^,  and 
Krupp  and  Cole?  showed  how  the  ideas  can  be  extended  to  unsteady  flow  with 
a  dimensionless  time  coordinate 

t  =  JJ  fi2/3t  (♦> 

c 

The  representative  point  (x,y,z)  runs  far  from  the  body  as  (6  ♦  0,  M„  ♦  1) 
for  fixed  (x,y,z)  to  express  the  fact  of  large  lateral  extent  of  disturbances 
near  M  *  1.  When  the  limit  process  expansion  is  substituted  into  the  basic  Euler 

OB 

m  m  • 

system  it  is  shown  that  to  this  order,  a  disturbance  potential  ♦(x,y,z,t)  exists 
such  that 

u  *  *x,  v  =  V*,  V  *  (!_.  »  i.)  (5) 

3y  3z 

This  disturbance  potential  vanishes  at  upstream  infinity  and  satisfies  the  basic 
TSD  equation 

(K  -  (y+1)$x)4>xx  +  “  2^xt  *  0  (6) 

The  pressure  coefficient  c  *  p~p«y .  is  found  from 

»  Tjm 


Cp  ■  -  2  62/3$x 


(7) 
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Some  properties  of  the  TSD  equation  are  now  summarized  for  steady  flow  in 
two-dimensions 

(K  -  (y+1)$x)$xx  +  ♦yy  a  0  TSO  eqn.  (8) 

The  equation  (8)  is  of  changing  type, 
elliptic  if  (y+1)4x  <  K 

hyperbolic  if  (y+1)4x  >  * 

with  flow  that  is  locally  subsonic  or  supersonic  respectively.  K  *  0  corresponds 

to  sonic  flow  at  infinity  (M  =  1  ).  The  local  Mach  number  M  can  be  shown  to  be 

•  l 

given  by 


The  typical  structure  of  flow  at  a  high  subsonic  Mach  number  is  shown  in  Fig.  3 
where  a  local  supersonic  (hyperbolic)  region  appears  over  an  airfoil.  The 
supersonic  region  contains  Mach  lines  or  hyperbolic  characteristics  given  by 


(y+1)*x-k 


(9) 


and  is  terminated  by  a  shock.  The  shock  conditions  are  contained  in  the 
conservation  form  corresponding  to  (8) 

(K*x  -  *0  (10) 

which  is  a  version  of  the  continuity  equation.  The  shock  is  a  discontinuity 
surface  across  which  $x,  $y  jump.  The  jump  conditions  are  given  by  the 

integrated  form  of  (10) 


in.) 


and  the  condition  that  4  is  continuous 
[♦]  *  0  or  Ctxldxs  ♦  [<fry]dys  ■  0 


(lib) 
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Here  (dxs,  dys)  are  line  elements  in  the  shock  surface. 

[  ]  *  jump  3  (  )b  -  (  )a  »  (  )b  *  quantity  behind  shock 

(  )a  =  quantity  ahead  of  shock 
Further  the  mass  flux  in  the  x-direction. 

pq  » 

*  »  1  +«*/3(k*  .Y+l  +  •••  *  (12) 

p«.U  *  2  * 

so  that  there  is  a  maximum  flux  at  local  sonic  speed  (c.f.  Fig.  4).  The  maximum 
of  the  mass  flux  corresponds  to  the  fact  that  stream  tubes  have  throats  at  local 
Mach  number  one.  The  TSD  equations  thus  contain  all  the  essential  features  of 
flow. 

III.  Early  GALC1T  Experiments 

A  small  transonic  wind  tunnel  2"  x  20*  was  constructed  in  1944-5  at  GALCIT 
by  Hans  U.  Liepmann  and  a  series  of  interesting  experiments  were  carried  out  in 
this  facility*.  Some  of  these  are  mentioned  here  and  some  In  later  sections. 
Experiments  were  carried  out  in  flow  past  a  series  of  circular  arc  airfoils  of 
dimensions  as  in  Fig.  5.  Surface  pressure  distributions  were  re<  „rded  and 
Schlieren  pictures  taken.  Fig.  6  shows  a  typical  pressure  distribution  at 
zero  angle  of  attack  (a  *  0).  The  local  supersonic  zone  is  shown  clearly  as  is 
the  substantial  difference  of  the  flow  with  laminar  and  turbulent  boundary  layers 
on  the  surface.  Turbulence  was  Induced  with  a  trip  wire.  For  the  turbulent  case 
In  which  the  boundary  layer  Is  thinner  the  shock  wave  terminating  the  sonic 
region  is  seen  clearly.  These  features  are  run  also  in  the  Schlieren  photographs 

t 

reproduced  in  Fig.  7. 

IV.  Numerical  Methods 

At  the  time  of  these  experiments  no  reliable  numerical  methods  existed  for 
calculating  the  ideal  flow,  although  Emmons  carried  out  some  relaxation 
calculations  for  flows  with  local  supersonic  regions.  Emmons'  method  however 
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was  inherently  unstable  and  did  not  resolve  shock  waves.  In  the  late  60's  Emmons' 
method  was  revised  to  eliminate  these  drawbacks*?  and  since  then  the  method  has 
undergone  substantial  development.  The  use  of  multigrids  has  speeded  convergence 
and  the  basic  ideas  have  been  extended  to  the  Euler  equations.  Here  we  describe 
the  original  idea  and  present  the  results  of  a  few  calculations.  Finite  dif¬ 
ference  methods  are  used  to  solve  numerically  the  boundary  value  problem  for  t 
corresponding  to  flow  past  an  airfoil.  The  best  results  are  obtained  by  using 
the  conservative  form  (8) 


(»*  -f'.Vraj =  o 


(13) 


and  a  corresponding  conservative  finite  difference  form.  The  essential  boundary 
conditions  are  (i)  tangent  flow  at  the  airfoil  surface 


(x.  0  ±)  *  FUjt  (x)  0  <  x  <  1 

where  y  *  6Fu>Jl(x)  represents  the  upper  and  lower  surfaces  respectively  (ii) 
vanishing  of  perturbations  at  infinity 


(14) 


4>x,  +  0  ^x^+y2 


(15) 


(iii)  Kutta  condition  that  the  flow  leaves  the  trailing  edge  smoothly.  In  TSD 
this  requirement  is  equivalent  to  zero  pressure  loading  at  the  tail  or  (c.f.7) 


♦x(1.0+)  =  ♦x(l,0-) 


(16) 


The  boundary  value  problem  is  shown  in  Fig.  8  where  in  addition  it  is  noted 
that  there  is  a  jump  in  $  across  the  wake 

O]  «  ♦(x,0+)  -  ♦(x,0-)  -  f  (17) 

where  r  =  ^)u  dx  +  v  dy  is  the  circulation.  It  is  also  noted  that  a«i  asymptotic 


-7- 


far  field  exists  which  for  subsonic  flow  has  the  form  of  a  circulation 

a  ♦  -  re  *  ...  ,  e  -  tan-1  (18) 

2iT  x 

Finite  difference  calculations  are  carried  out  on  the  (i»j)  mesh  indicated.  A 

conservative  form  is  derived  by  considering  fluxes  around  a  central  box.  This 

reasoning  is  extended  so  that  the  shock  waves  which  appear  spread  over  three  or 

four  mesh  points  are  consistent.  For  stability  the  difference  scheme  must  be 

chosen  to  be  type  sensitive.  The  solutions  at  (i,j)  can  be  influenced  only  by 

upstream  points  if  the  flow  is  locally  supersonic  ($x  >  K/(y+1))»  but  by  both 

upstream  and  downstream  points  if  the  flow  is  locally  subsonic.  In  the  finite 

difference  approximations  $x  at  (i,j)  can  be  calculated  from  a  centered  formula 

$(c)  involving  (i+1,  1-1 )  or  a  backward  formula  $(&)  involving  (i,i-2).  When 
x  x 

these  agree  an  (i,j)  can  be  designated  subsonic  or  supersonic  as  indicated  in  Fig. 
9.  At  subsonic  points  (13}  provides  an  explicit  equation  for  in  terms  of 
neighbors  on  all  sides,  the  computational  star  in  Fig.  9.  For  supersonic  points 
however  an  implicit  scheme  is  used  involving  only  upstream  points.  Two  other 
kinds  of  points  denoted  as  "sonic",  which  applies  for  points  near  the  sonic  line 
where  the  flow  accelerates  and  "shock"  which  applies  for  points  where  the  flow 
decelerates  to  subsonic  through  a  shock  are  shown. 

Since  the  problem  is  non-linear  the  local  state  is  not  known  in  advance  and 
an  iteration  scheme  is  used.  At  a  given  iteration  the  difference  schemes  are 
chosen  and  the  unknowns  solved  for  on  a  vertical  line.  Sweeps  are  made  in  the 
downstream  direction.  The  analytical  far  field  is  used  as  a  boundary  condition  but 
the  value  of  F  is  adjusted  as  the  trailing  edge  is  passed.  To  speed  convergence 
tf|e  latest  values  are  used  whenever  possible.  The  method  converges  well,  gives  the 
correct  shock  jumps,  and  for  lifting  cases  automatically  satisfies  the  Kutta 
condition. 
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Some  results  of  calculations  are  shown  in  Fig.  10  and  11  for  a  parabolic 
arc  airfoil 9  and  a  NACA  001210.  For  the  first  the  sonic  line  and  shock  appear 
as  well  as  the  surface  pressure.  Figure  12^®]  shows  the  flow  field  features  for 
a  higher  freestream  Mach  number  where  the  main  shock  has  moved  aft  of  the  airfoil 
and  a  fishtail  shock  pattern  appears. 

V.  Shock-free  Supercritical  Flows 

The  possibility  of  shock  free  flows  past  airfoils  with  local  supersonic 
regions  was  of  considerable  interest  in  the  early  days  of  transonic  research. 

This  was  stimulated  by  Ringleb's  exact  hodograph  solution*!,  analogous  to 
incompressible  flow  around  a  half-plane,  which  has  a  smooth  transonic  region. 
Experimentally  it  was  possible  to  produce  a  small  shock-free  supersonic  zone  in 
the  flow  around  a  simple  shape.  For  example  see  Fig.  13  where  some  of  Liepmann's 
results  are  reproduced4. 

Other  exact  solutions,  analytical  and  numerical,  were  derived  from  hodograph 
considerations  and  gave  shock  free  flow  past  special  airfoil  shapes.  The  hodograph 
equations  are  linear  so  that  solutions  can  be  obtained  and  the  airfoil  shape  foun 
later.  For  the  TSD  system  the  hodograph  equations  are  obtained  by  a  direct 
interchange  of  dependent  and  independent  variables.  Rewrite  (8)  as  the  system 


where 

w  *  (y  +1  )$x  -  K  ,  v  =  (y+1  )4»y 


(19) 


Then 


wyv  =  xw 
xv  *  yw 


(20) 


since  wx 


etc,  and  J 


Jacobian  «  .  From  this,  it  is  seen 

3(w,v) 
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that  the  approximate  stream  function  y(w,v)  is  given  by  a  solution  of  Tricomi's 
equation. 

wyw  “  Pww  =  0  (21) 

This  is  the  simplest  linear  equation  of  mixed  type,  elliptic  in  the  subsonic 
region  w  <  0,  hyperbolic  in  the  supersonic  region  w  >  0.  The  exact  potential 
equation  can  be  transformed  in  a  similar  way  to  produce  Chaplygin's  equation  which 
has  properties  analogous  to  (21). 

Some  special  exact  solutions  for  mixed  flow  past  an  airfoil  were  given  by 
Tomotika  and  Tamadal2.  Recent  advances  have  been  the  work  of  Nieuwland1^ 
which  used  Chaplygin  functions  to  represent  analytically  families  of  airfoils  and 
Garabedian  and  Korn14  who  use  a  finite-difference  hodograph  method  to  obtain 
numerical  solutions  of  airfoil  flows.  The  latter  employ  real  characteristics  in 
the  hyperbolic  region  and  complex  characateristics  in  the  elliptic  region  and  a 
sufficient  number  of  parameters  to  generate  families  of  shapes  with  local  subsonic 
zones.  A  typical  airfoil  and  its  pressure  distribution  in  shock  free  flow  is 
shown  in  Fig.  14.  The  pressure  distribution  is  also  calculated  according  to  TSD 
theory  with  a  fully  conservative  relaxation  scheme  ( FCR)  as  described  in  Sec.  IV 
and  a  non-conservative  scheme  (NCR)16,  as  well  as  a  calculation  of  Garabedian, 

Korn  &  Jameson16.  The  drag  coefficient,  which  is  theoretically  zero  (see  next 
section)  is  evaluated  on  several  different  control  surfaces16. 

The  smooth  hodograph  solutions  are  seen  to  be  isolated  solutions  since  a 
smooth  mixed  flow  can  not  be  found  when  the  boundary  conditions  in  the  physical 
plane  change  slightly  (Morawetz1?).  An  early  paper  of  Guderley16  suggests  a 
singularity  of  a  perturbation  in  the  downstream  corner  of  the  sonic  region.  They 
are  isolated  solutions  in  the  same  sense  as  the  famous  Busemann  supersonic 
biplane.  When  the  conditions  are  changed  slightly  a  neighboring  solution  is  found 
with  a  shock  wave.  An  example  of  the  calculation  of  such  a  flow  appears  in 


P 
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Fig.  15,  where  the  shock  wave  is  apparent.  Experiments  verify  these  features. 

Guderley^  suggests  that  the  shocks  in  the  local  supersonic  region  are 
initiated  by  an  envelope  formation  of  the  compressive  reflection  of  Mach  lines 
from  the  sonic  line.  A  detailed  look  at  the  calculations  of  Fig.  10  verify  this 
idea^.  Thus  a  shock  free  airfoil  is  one  in  which  the  shape  of  the  sonic  line 
is  just  right  so  that  no  envelope  is  formed. 

VI.  Shock  Waves  and  Drag 

The  connection  of  shock  waves  and  drag  was  made  explicit  by  Liepmannl^ 
in  1950  with  a  study  of  linearized  supersonic  flow  past  an  airfoil.  By 
considering  the  weak  shock  waves  in  such  a  flow  a  correction  to  the  Mach  angle  was 
made  for  shock  angle  and  location  and  the  formula  was  found 

Drag/ (length)  =  p«T„/[S]dy  (22) 

shocks 

where  [S]  *  jump  in  specific  entropy  across  a  shock  wave.  Drag  is  related  directly 
to  entropy  production.  Similar  results  appear  later  in  works  by  von  Karman^O 
and  0swatitsch21. 

Analogous  considerations  apply  to  TSD  flow.  Germain22  gave  a  derivation  of 
a  drag  formula  for  TSD  theory  starting  from  a  local  conservation  law  in  two- 
dimensions.  The  generalization  to  three-dimensional  flow  was  given  in5  and 
reads 


r 


This  conservation  law  easily  follows  from  the  three-dimensional  version  of  (8) 
and  the  equation  of  irrotational ity.  The  values  of  (uv)  on  y  *  0  are 

proportional  to  the  incremental  drag  for  a  planar  system  since  the  pressure 

increment  is  proportional  to  u  and  the  wing  slope  to  v.  Integration  of  (23)  over 


/ 
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a  control  surface  as  shown  in  Fig.  16  yields  a  formula  for  drag.  Since  (23)  is 
not  a  physical  global  conservation  law  it  is  not  conserved  across  shock  waves. 
Letting  the  control  surface  grow  to  infinity  contributions  from  the  shock  and 
vortex  sheet  remain.  The  result  is 

D  = _ =  lim  //  (v£+w£)dy  dz  -  Ill  //  [u]3dy  d£  (24) 

p  u2<S4/3c2  x‘*’e#  2  -12 

*  vortex  drag  shocks 

The  first  term  of  (24)  is  the  familiar  expression  for  vortex  drag  in  terms  of  the 

kinetic  energy  in  the  wake  while  the  second  term,  the  wave  drag  is  a  scaled 

version  of  the  Liepmann  formula,  [u]  is  the  jump  of  $x  across  the  shock  and  the 

entropy  change  is  proportional  to  [u]3.  Thermodynamic  considerations  for  weak 

shocks  show. 

II  [u]3dy  dz  ~  p.T.  //  [S]dy  dz  (25) 

shocks  shocks 

It  should  be  noted  that  these  simple  results  do  not  carry  over  for  stronger 

shocks.  The  drag-entropy  formula  (24)  can  be  used  to  provide  a  check  on  the 

consistency  of  numerical  TSD  calculations^.  jn  Fig.  (17)  the  shock  and  flow 

field  for  a  parabolic  arc  at  H,  =  .909  is  shown.  The  J[u]3dy  -  /[cp]3dy  is 

plotted  also  and  the  drag  coefficient  is  computed  both  from  the  surface 

pressure  integral  and  the  entropy  jump.  The  result  is 


Cd  =  .0315  surface  pressure  integration 
Cd  =  .0320  shock  entropy  integration 


(26) 
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VII.  Sonic  Flow  and  the  Law  of  Stabilization 

The  flow  at  M»  =  1  has  a  special  structure  elucidated  first  by  Guderley  and 
Frankl .  In  order  to  understand  this  flow  consider  first  the  sequence  of  flow 
patterns  past  an  airfoil  or  body  at  free  stream  Mach  numbers  close  to  one.  Fig. 

19.  At  high  subsonic  Mach  numbers  a  large  supersonic  zone  is  formed  terminated 
by  an  oblique  shock  near  the  trailing  edge  and  the  fishtail  shock,  as  calculated 
in  Fig.  12.  At  M»  =  1,  the  supersonic  zone  in  the  sequence  of  steady  flows  grows 
to  infinity  and  the  fishtail  shock  goes  to  downstream  infinity.  A  limit  charac¬ 
teristic  or  Mach  wave  appears  which  is  asymptotically  parallel  to  both  the  sonic 
line  and  the  tail  shock  at  infinity.  This  limit  characteristic  divides  the  flow 
field  into  an  upstream  and  downstream  part.  Any  (infinitesimal)  disturbance  in 
the  flow,  in  the  supersonic  region,  for  example,  can  send  a  disturbance  downstream 
which  eventually  reaches  the  sonic  line  and  affects  the  entire  subsonic  region. 

Any  disturbance  originating  downstream  of  the  limit  characteristic  can  not  affect 
the  upstream  flow.  The  flow  in  the  upstream  section  up  to  the  limit  character¬ 
istic  thus  has  an  elliptic  nature  and  must  be  calculated  all  at  once.  It  is 
effectively  independent  of  the  flow  downstream  of  this  characteristic.  Downs  ;rsam 
of  the  limit  characteristic  the  flow  can,  for  example,  be  calculated  as  in  a 
hyperbolic  region  by  the  method  of  characteristics.  When  the  free  stream  becr-mss 

slightly  supersonic  a  detached  shock  wave  which  has  subsonic  flow  behind  nerr  the 
axis  (y=0)  appears  ahead  of  the  body.  The  flow  becomes  supersonic  near  the 

airfoil  and  terminates  again  in  an  oblique  tail  shock.  It  can  thus  be  noted  that 

the  flow  in  the  neighborhood  of  the  body  does  not  change  qualitatively  very  much 

t 

since  the  oncoming  flow  is  always  close  to  a  uniform  subsonic  flow. 

For  M»  «  1,  the  far-field  is  a  similarity  solution.  With  K  «  0  (8)  becomes, 
for  the  planar  case 


(y+1)$x$xx  "  $yy  a  0 


(27) 
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The  far-field  can  be  thought  of  as  being  produced  by  a  singularity  at  the 
origin  and,  due  to  non-linearity,  represented  in  the  form 

(ir+l)4>(x,y)  *  y3<‘2  f(5;<),  5*4 

y* 

a  one-parameter  family  of  similarity  solutions.  '  Thus  note 

(Y+l)*x  =  p2<-2  f(5)  ,  (y+l)$xx  =  y<‘2  f-(C)  (29) 

f  ( 5  ;< )  is  thus  found  from 

(f '-ie2$2)f"  .  5<(l-tc)5f '  +  3(l-<)(3<-2)f=0,  -  -  <  5  <  5l  (30) 

The  coordinates  are  illustrated  in  Fig.  20.  The  flow  must  start  at  infinity 
decelerating  and  spreading  around  the  body,  then  accelerate  through  sonic  and 
smoothly  through  the  limit  characteristic.  It  can  be  seen  from  (29)  that  the 
sonic  line  (f'=0)  as  well  as  the  limit  characteristic  and  tail  shock  must 
lie  on  similarity  curves  5,5l»  and  5s  respectively.  The  characteristic 
condition  (+)  (29)  becomes 


f '  (cl)  * 


(31) 


4  also  must  be  symmetric  in  y,  the  apparent  thickness  ev”sc:  dominatrr  toe 
flow  at  infinity.  This  occurs  because  both  acceleration  and  decelerr.  .io:i  of  a 
sonic  flow  produces  a  widening  of  stream  tubes  (cf  Fig.  4),  One  con: tar.-; 
of  integration  remains  if  integration  is  started  at  5  *  -  •  and  two  conditions 

t 

must  be  satisfied  at  the  limit  characteristic  if  f"(5L)  is  to  be  finite.  Thus 
the  solution  exists  only  for  a  special  of  k.  Guderley  snowed  from  trs 
hodograph  solution  that 


K 


4 

I 


planar  case 


(32) 
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and  by  numerical  Integration  that  k  *  A  for  the  analogous  axisymmetric  case. 
The  latter  result  was  shown  to  be  exact  by  various  writers  e.g.23.  Eqn.  (30) 


has  a  group  property  so  that  the  solution  can  be  written 


f(C)  =  1_  F(AC),  A  =  const.  (33) 

A3 

The  scale  factor  A  depends  on  the  size  (and  shape)  of  the  body.  Thus  the 
far-field  can  be  standardized  to  have  Cl  -  1*  By  considering  shock  jumps  the 
solution  can  be  continued  to  downstream  infinity.  A  particularly  useful 
representation  of  the  solution  was  given  parametrical ly  by  Frankl24 


n.) 

1  55 

C(s)  -  ai1/5s-2/5  -  lj 

where  the  following  relationships  are  noted 

5  0  C*  CL  *  1 

s  Oil  A/3 

2 


Cs 


1  (5/1  +  8) 
6 


(35) 


Frankl  produced  this  result  by  clever  observations  about  special  hodograph 
solutions;  an  analogous  result  has  been  derived  by  inspection  in  the 
axisymmetric  case. 

A  useful  extension  of  the  far-field  (28)  is  to  regard  it  as  the  first  term 
of  an  expansion  of  the  form  , 


(Y+l )*(x,y)  *  y2/5  i-  f(AC)  +  C  y°°  Mi  +  C  yal  liiMl  (36> 

A3  u  A3  1  A3 

valid  as  y  ♦  •  for  C  fixed.  The  perturbation  functions  fj  >  g,  j  >  0.  1 
satisfy  the  variational  equation  (with  *  o) 
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(f '  -  16  C2)g"  +  (fH  +  A  (2a  -  I)$)g'  -  a(a  -  l)g  «  0  (37) 

25  5  5 

If  g  =  s~a/2  h(t),  t  *  3  s  (37)  becomes  a  standard  hypergeometric  equation 

T 

t(l-t)hH  +  (I  -  i  (5a  ♦  4)t)h *  +  1_  a(5a  -  2)h  -  0  (38) 

2  3  12 


in  the  interval 

around  t  =  0  are 


\  0  <  t  <  I  / 


The  two  linearly  independent  solutions 


hj  -  F(5  a,  1  (2-5a) ;  l;t) 


l 

hH  =  t2  F(1  (5a+l),  |  (1-a);  |;t) 


(39) 


Thus  solutions  smooth  on  passing  through  the  limit  characteristic  are  given  by 
the  spectrum 


o . 3,  _  2  ,-1,0,2  8  11  14 

i  *  r  i  5  ’  ’  5  ’  t  ’  r  *  r 


(40) 


For  the  expansion  near  infinity  a  -  0,  o  =  -  1  and  tte  f*  are  simole 

0  1  ■5- 

functions. 

The  solutions  just  discussed  can  be  used  to  relate  a  flow  with  a 
free-stream  Mach  number  close  to  one  to  that  at  Mach  number  one.  In  Tf?  form 
this  relates  a  flow  with  small  |K|  to  that  with  K  *  0.  For  flow  subsor at 
infinity  the  dominant  term  in  the  far-field  is  the  circulation  term,  basically 
the  solution  which  decays  most  slowly  for" 


K  A  +  ♦ 
*xx  *yy 


0  ,  ♦  ♦  -  i-  e 


(41) 


But  for  K  ■  0  the  far-field  is  given  by  (26)  with  k  *  1  .  Thus  there  is  a  non- 

5 

uniformity  at  Infinity  of  ♦(x,y;K)  as  K  +  0  and  inner  and  outer  expansions  are 
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needed.  For  the  Inner  expansion  valid  near  the  airfoil  the  first  term  is  sonic 
flow  past  the  airfoil  and  corrections  are  sought.  In  early  work  Liepmann  and 
Bryson28  proposed  that  at  M»  *  1  the  local  Mach  number  on  a  body  does  not 
change  as  the  free  stream  Mach  number  changes  near  one  because  of  the 
qualitative  ideas  outlined  above.  The  Law  of  Stabilization  proposed  by 
Ryzhov  and  Lipschitz28  and  the  more  detailed  work  of  Cook  and  Ziegler2?  following 
the  method  of  matched  asymptotic  expansion  indicated  in28  give  a  deeper  and 
more  precise  description  of  this  Mach  number  freeze.  The  starting  point  is  a 
solution  at  M»  *  1  and  several  are  available,  for  example,  for  a  wedge  Guderley 
and  Yoshihara8  gave  the  result  while  Tse28  worked  out  lifting  airfoils  at 

M»  *  1.  Both  of  these  methods  rely  on  a  hodograph  formulation.  Let  ♦*(x,y) 

represent  the  sonic  flow  past  an  airfoil  with  the  far-field  (28).  Then  the  flow 
for  K  >  0  is  represented  as  an  inner  expansion,  generated  from  a  limit  process 

K  -*■  0,  x,y  fixed. 

*(x,y;K)  -  **(x,y)  +  K  *  +  e(K)  ♦  (x,y)  +  ,M  (42) 

Y+l  t 

e(K)  is  the  order  of  magnitude  of  the  correction  which  is  sought  and  ♦c(x,y) 

is  the  correction  potential.  $*(x,y)  satisifes  the  sonic  TSD  eqn  (25)  and  its 

boundary  conditions  while  $c(x,y)  is  found  from  the  variational  equation 

(y+1)U*  4»c  +  ♦*  4c  )  -  ♦c’“  *  0  (43) 

x  xx  xx  x'  yy 

e 

The  expansion  (42)  is  not  valid  near  infinity  so  that  an  outer  expansion  in 
rescaled  variables  is  sought 

♦  (x,y;K)  *  o(K)  4(x,y)  +  (44) 
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The  associated  limit  process  has  x,y  fixed  as  K  +  0,  x,y  ♦  •  ,  that  is 
x  =  u(K)x,  y  *  v(K)y  p,  v  -*•  0. 

The  condition  that  the  rescaled  equation  be  subsonic  at  infinity  but  still 
nonlinear  and  of  changing  type  produces  a  one  parameter  family  of  flows. 

♦(x,y;K)  «  K.  $(x,y)  +  with  v  *  p/K‘  (45) 

p(K) 

with  the  resulting  version  of  the  TSD  equation 

(1  -  (y+1)$x)$xx  +  fyy  *  0  (46) 

The  inner  and  outer  expansions  must  match  and  in  a  simple  way  as  x,y  ♦  •  , 
x,y  >  0.  is  defined  along  similarity  curves  which  then  must  match 

£  =  x  s  K2/5  x  x  r  . 

{  *4/5  ^75-  ;  4/5  « 

Thus  p  =  K2,  x  ■  K2  x,y  *  K5/2  y  .  The  general  form  of  the  flow  as  it  looks  in 
outer  variables  is  shown  in  Fig.  21.  Higher  terms  in  the  flow  near  infinity 
for  the  inner  and  the  flow  near  the  origin  in  the  outer  expansion  (42), (45) 
then  depend  on  the  spectrum  associated  with  (37). 

These  expansions  can  then  be  written  (with  the  scale  factor  A  *  1) 

♦  (x,y;K)  «  _I_|y2/5  f(c)  +  cQ  +  }‘1/5  fjU)  +  ’ ' ' J  ♦  J&j. 

+  e(K)*c(x,y)+  <47a> 

inner 

»  |  _l_jy  2/5  fl5>  +  *  +  ci  y  8/5/Ti(|)  +  ’“j,  outer  (47b) 

A  comparison  of  these  two  as  x,y  ♦  0,  x,y  ♦  •  along  similarity  curves  shows  that 
♦l(x,y)  has  similarity  form  at  infinity  and  that 
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e(K)  =  K3 


(48) 


This  weak  dependence  on  deviation  from  H»  =  1,  for  the  potential  and 

pressure  distribution  is  the  essential  part  of  the  law  of  stabilization.  The 
term  k  x  in  these  expansions  represents  the  change  of  the  flow  at  infinity 
and  adds  a  constant  pressure  level  to  the  solution.  Germain22  derived  a 
conservation  law  formula  for  the  scale  factor  A  relating  it  to  properties  of  the 
solution  (not  known  in  advance)  on  the  body  surface.  Cook  and  Ziegler27  have 
extended  these  ideas  to  find  Cj.  The  boundary  value  problem  for  the  inner 
correction  solution  is  sketched  in  Fig.  22.  An  empirical  fit  of  the  Law  of 
Stabilization  to  some  experimental  drag  measurements  of  Vincenti3®  is  drawn  in 
Fig. 22. 


VIII.  Concluding  Remarks 

TSD  theory  ard  calculations  have  made  great  studies  in  the  last  35  years 

based  on  the  work  of  the  pioneers.  The  understanding  of  the  physical  phenomena 

that  has  been  thus  achieved  enables  us  to  go  forward  with  more  elaborate 

calculations  and  including  more  physical  effects.  One  can  only  admire  the  deep 

insight  of  the  early  workers  in  this  area  such  as  H.W.  Liepmann. 
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